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Abstract 

A recursive calculational scheme is developed for matrix elements in the generalized 
seniority scheme for the nuclear shell model. Recurrence relations are derived which 
permit straightforward and efficient computation of matrix elements of one-body and 
two-body operators and basis state overlaps. 
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1. Introduction 

The generalized seniority [1, 2] or broken pair [3, 4] approximation provides a trunca- 
tion scheme for the nuclear shell model, based on the dominance of like-nucleon pairing 
effects in semimagic (or nearly semimagic) nuclei. The fundamental premise is that there 
exists an energetically favored collective S pair, which can be constructed from like nu- 
cleons coupled to angular momentum zero, and that the ground state of an even-even 
nucleus can be well approximated by a condensate of such collective pairs. Low-lying ex- 
cited states may then be obtained by breaking one or more S pairs. The resulting states 
are classified by the generalized seniority v, defined as the number of valence nucleons 
not participating in a collective S pair. The generalized seniority approach effectively 
reduces the dimensionality of the n-particle valence shell model to that of a w-particle 
shell model problem. In practice, such generalized seniority shell model calculations have 
been carried out for v < 4 (e.g., Refs. [5-9]). 

Moreover, the generalized seniority scheme provides a microscopic foundation for the 
phenomenologically successful interacting boson model (IBM) [10] and, for odd-mass and 
odd-odd nuclei, the interacting boson fermion model (IBFM) [11]. In this context, gen- 
eralized seniority states constructed from the S pair and a collective D pair (consisting 
of like nucleons coupled to angular momentum 2) are mapped onto IBM states built 
from analogous combinations of s and d bosons [12-16]. The need for a fully microscopic 
derivation of the IBM Hamiltonian parameters and transition operators is a longstand- 
ing problem. A well-developed microscopic mapping would be especially valuable for the 
IBFM, where the relevance of single-particle degrees of freedom is particularly manifest. 
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For instance, higher-order mapping from the generalized seniority shell model space onto 
the IBFM would facilitate the prediction of collective effects in /3 decay [17]. The need 
is also emphasized by recent applications of the model to neutrinoless double beta de- 
cay [18], where phcnomenological calibration of the model is not a viable alternative. 
For microscopic derivation of the IBM and IBFM, matrix elements of various operators 
between states of generalized seniority v < 6 are of interest. 

In addition to matrix elements of operators between generalized seniority states, it 
is also necessary to compute the overlaps of these states, since the canonical construc- 
tion [3] leads to a nonorthogonal basis. The S pair creation operator is obtained as a 
linear combination of pair creation operators for each active single-particle level (i.e., j 
shell). The individual j shells may be treated using the conventional seniority or qua- 
sispin results [19, 20]. Then, the overlaps and matrix elements in the generalized seniority 
scheme may in principle be computed from the single-shell results by combinatorial ar- 
guments [3, 13, 21]. 

However, the combinatorial derivation [13] rapidly becomes cumbersome with in- 
creasing generalized seniority. It was noted [22] that commutator methods can be used 
to obtain simplifications, but the resulting hybrid approach could still only be practi- 
cally applied for v < 2 [22, 23]. The main challenge in proceeding to higher generalized 
seniority lies in accounting for the intermediate angular momentum couplings which oc- 
cur among the particles not participating in the collective S pair. Several approaches 
have been pursued. The combinatorial approach may be pushed further by introducing 
certain simpler intermediate quantities, essentially norms for uncoupled states, and then 
reexpressing the overlaps for the angular-momentum coupled states in terms of these 
(however, matrix elements were not explicitly considered in this formulation) [24]. One 
may directly carry out the calculations of overlaps and matrix elements by Wick's the- 
orem in an uncoupled scheme, subsequently recoupling the results [25]. An alternative, 
indirect construction of the generalized seniority basis relies upon number projection of 
generalized-quasiparticle states involving complex parameters. Matrix elements are then 
evaluated as contour integrals in the complex plane [4, 5, 26]. These latter methods have 
been applied for v < 4. 

In the present work, a systematic calculational scheme is established, in which recur- 
rence relations for the matrix elements and overlaps are derived by angular-momentum 
coupled commutator methods [27, 28]. This is essentially a Wick's theorem approach, 
which, however, retains the angular-momentum coupled structure. General expressions 
are obtained by which the one-body operator matrix elements and the overlaps for states 
of generalized seniority v and pair number TV are expressed in terms of commutators 
which yield matrix elements and overlaps involving lower values of v and/or N. The 
explicit recurrence relations to be used in the numerical calculations for a specific value 
of the generalized seniority are obtained from these generic expressions by a straight- 
forward procedure involving repeated application of a commutator product rule. The 
matrix elements of two-body operators can be computed directly in terms of these basic 
one-body operator matrix elements and overlaps. An essential motivation for the present 
approach is that calculation can proceed to higher generalized seniority in a methodical 
fashion, through systematic application of this coupled commutator product rule, and 
that the process is readily amenable to automation [29]. For illustration, the explicit 
recurrence relations are given for matrix elements and overlaps involving low generalized 
seniority states, sufficient for generalized seniority shell model calculations with v < 3. In 
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particular, calculation of the matrix elements of a two-body Hamiltonian between states 
of v = 3 requires one-body operator matrix elements for v < 4 and overlaps for v < 5. 

After a review of definitions (Sec. 2), the general commutation scheme is outlined 
(Sec. 3), and the underlying commutator algebra results are summarized (Sec. 4). A 
coupled set of recurrence relations for the one-body operator matrix elements (Sec. 5) 
and overlaps (Sec. 6) are obtained. Matrix elements of two-body operators are then 
expressed in terms of these quantities (Sec. 7). 

2. Definitions 

First, let us review the definitions necessary for the generalized seniority scheme. 
Within a single major shell, the angular-momentum label is sufficient to uniquely specify a 
single-particle level, so let C\ 7 be the creation operator for a particle in the state of angu- 
lar momentum c and ^-projection quantum number 7. 1 The angular-momentum coupled 
product of two spherical tensor operators is defined by (A a xB b ) c i y = ^Z a p{aoibj3\c^)A'^B b p, 
and we follow the time reversal phase convention A a a = {— ) a ~ a A a _ a [30]. Then the 
angular-momentum coupled pair creation operator is defined by 

Kl = (ci x c\r, (i) 

and its time-reversed adjoint is A e ab = (A^l) = —(C a x Cb) e ■ 

The state of zero generalized seniority is defined by the S'-pair condensate, \S N ) = 
S'^ Ar |0), where the collective S pair in the generalized seniority scheme is defined by 

St = £a c |^ c t, (2) 

c 

with c = (2c + f ) 1 / 2 . The amplitudes a c are conventionally taken subject to the nor- 
malization condition J2 C (^ C + ^) a c = J2 C (^ C + -0- They may be obtained by a variety 
of prescriptions [1, 3, 25, 31, 32], for instance, variationally so as to minimize the energy 
expectation value (S N \H\S N ). 

More generally, a state of generalized seniority v is constructed as \S N F :f ) — 
N F^\0), where the fermion cluster is a product of v creation operators, cou- 
pled to total angular momentum /. A complete set of such clusters is constructed 
by successive couplings of the form = {{A^ x A^) ■ ■ ■ x ^4^)^ for v even or 
(((A'f x A^) ■ ■ ■ x A^) x C^y for v odd. In particular, the conventional microscopic inter- 
pretation [14] of the IBM is formulated in terms of generalized seniority states involving 
the collective D pair, defined as a general linear combination 

ab 
a<b 



1 For simplicity of notation, we do not distinguish between the level c and its angular momentum j c . 
Also note that, following Ref. [28], we denote operators by capital roman letters, angular momenta by 
lower case roman letters, and angular-momentum ^-projection quantum numbers, when needed, by the 
corresponding lower case Greek letter (e.g, A%, B^, C°). 
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of pairs of angular momentum 2. The states of interest for mapping from the shell model 
to the IBM, through the Otsuka-Arima-Iachello (OAI) mapping [12-16], are of the form 
\S N ), \S N D), \S N (DD)f), etc. For the analogous mapping [33] to the IBFM, the relevant 
states are \S N C a ), \S N (DC a y), etc. Any matrix element or overlap involving D pairs 
may be expanded in terms of those for the elementary pairs A 2 ab , through (3). 

The matrix elements of most immediate interest are those of one-body and two-body 
operators. Any spherical tensor one-body operator may be expressed in terms of the 
elementary one-body multipole operators 

Tl s = [C\ x C s )\ (4) 

through the usual second-quantized realization, which in angular-momentum coupled 
form becomes U u = -^^Z^u" 1 {a\\U u \\b)T^ b . 2 Thus, we must consider matrix elements 
of T' s taken between states of the generalized seniority scheme, which are of the form 
(S N G 9 \\T* S \\S N F^), where and G 9 represent clusters of v nucleons not participating 
in the collective S pairs. Since the generalized seniority states are not orthonormal, as 
canonically constructed in terms of clusters above [15], it is also necessary to compute 
the norms and overlaps (S N G f \S N Ff). 

The creation operators for the generalized seniority basis states, as defined above, 
have a uniform structure consisting of sequentially coupled pair creation operators 
(with a single additional creation operator in the case of odd v). Therefore, it is 
more compact and readable to simply label these states by the single-particle level la- 
bels and the angular momenta for intermediate couplings, that is, as \S N ) = S^ N \Q), 
\S N c) = & N C%\0), \S N {ab) e ) = S^ N A e J\0), \S N (ab)H 9 ) = N (A e J x Cj)»|0>, 

^(abY^r 9 ) = & N {A e J x A%i)9\0), etc. 

Only matrix elements and overlaps involving bra and ket states of equal generalized 
seniority need be calculated explicitly through recurrence relations. To facilitate writing 



uuuee nie symuui (_ 

for the overlaps of states of equal generalized seniority v, as 



(and discussing) the recurrence relations, we furthermore introduce the symbol O n 



= (S N \S N ) 

lW^U = /qNiqN, 



0%>[c\c] = {S»c\S»c} 
l) e \{ab) e ] = (S N {cd) e \S 
0%>[(cdy f\{ab)H°] = {S N {cdy ] 9 \S N {abfi 9 ) 



0^[{cdy\{ab) e ] = (S N (cd) e \S N (aby) (5) 



0f[{cd) f {kl) n9 \{ab) e {ij) m9 ] ee (S N (cdy\kl) n9 \S N (ab) e (ij) m9 ), 

etc. Similarly, we introduce T^[- ■ ■} for the one-body operator matrix elements between 
these states, e.g., 

T^ ) [( C( i) / (fc/) rl ' l |(r S )*|(a6) e (ij) m 9] ee (S N (cdy(kl) nh \\T* s \\S N {aby(ij) m9 ). (6) 

These definitions are intended to follow the bracket notation as closely as possible, while 
clearly exhibiting the labels v and N with respect to which the recurrence will be carried 



2 Howcver, in comparing with, e.g., Ref. [34], note that the overall sign of this expression varies 
depending upon the time-reversal phase convention, presently A™ = (— ) a ~ a A a ^ a . 
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out, and also more simply laying out the arguments of and T^' considered as 
symbols in the recursive computational scheme. 

Matrix elements or overlaps involving states of unequal generalized seniority can read- 
ily be expressed in terms of these by expanding one or more S pairs in terms of the pairs 
(aa)°, i.e., created by A° a \, for individual j shells, according to the definition (2), e.g., 

<s"(cd)»is"-v) e (yr fl > = £«^ ( ? ) 

k 

Here we notice that, although the pair (kk)° carries seniority zero, it is not the collective 
S pair and therefore carries a generalized seniority of 2. 

It should be noted that the and 0$ obey a variety of symmetry relations 

under rearrangement of the arguments, and therefore they need not all be calculated 
independently. For instance, different ordcrings of single-particle labels within the bra or 
ket are related by, e.g., A\l = -6(abe)A e a l and (A e J x A™^) 9 = e(emg)(A^ x A e Jy, 
where 

e(abc---) = (-) a+b+c+ -. (8) 

Thus, for instance, 

T%\cdy{kl) n \rsy\{aby{ij) m °] = -e(abe)T^\cdy{kl) n \rs)%baf(ijr% (9) 
interchanging particles within a pair, as indicated by the arrows, or 

T^ 4) [(cd)/(fcO n ' l |(^) t |(a6) e fe-) m9 ] = 6{emg)T%\cdy{kiy%rsY\(ijT{aby9], (10) 

< > 

interchanging pairs within a cluster. Moreover, the property of the reduced matrix 
element under complex conjugation, (c||£? fc ||a)* = (— ) a+fc ~ c (a||iT' i '||c), gives the relation 

(S N G g \\T* s \\S N F f ) = -(-y+*+9-f(s N F f \\T t sr \\S N G g ), (11) 
where we note that C\ = —C\, and consequently T*t = —0(rst)T* r . Thus, for instance, 
T^\cdy{kiy\rs)%aby{ijr^] = -6{rsgh)T^ 

i > 

(12) 

interchanging bra and ket and conjugating the operator. 

3. General scheme for recurrence 

Recurrence relations for the reduced matrix elements and overlaps can be naturally 
derived by commutator methods. These matrix elements and overlaps must first be ex- 
pressed as vacuum expectation values of coupled products of creation and annihilation 
operators. The operators are then reordered, by making use of their commutation (or 
anticommutation) relations, until all the resulting terms themselves represent matrix 
elements or overlaps of generalized seniority states. Since the commutators (or anticom- 
mutators) yield terms with fewer total creation and annihilation operators, the matrix 
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elements and overlaps resulting from these terms will automatically involve states of lower 
pair number or generalized seniority, thus giving rise to relations which are recursive with 
respect to JV or w. 

In general, consider two states \A^) — A^\0) and |C^) = Cf;t|0), obtained by the 
action of some operators A^J and on the vacuum, for instance, coupled products of 
fcrmion creation operators. Then the reduced matrix element of an operator B b can be 
reexpressed as the vacuum expectation value 

(C c \\B b \\A a ) = (-) a - b - c (0|[C' c x (B b x A at ) c ]» (13) 

of a scalar triple product of operators. 3 The overlap (B a \A a ) = (B^\A^,) of two states 
can similarly be evaluated as the vacuum expectation value 

(B a \A a ) = a- 1 (0\(B a x A at )». (14) 

The matrix element of the one-body operator T* s between states of generalized se- 
niority v is therefore given by an expression of the form 

(S N G 9 \\T t rs \\S N F f ) = (-)*-*-<> {0\(G a S N x T l rs x S* N F f *)°\0), (15) 

where t and G 9 t represent two clusters consisting of v fermionic creation operators. 
The overlap of two states of generalized seniority v is similarly 

(S N G f \S N F f ) = f- 1 (0\{G f S N x & N F f *f\0), (16) 

where again F* t and t represent two clusters consisting of v fermionic creation oper- 
ators, now with the same angular momentum / = g. 

To outline the general scheme for deriving the recurrence relation for a matrix ele- 
ment , let us momentarily suppress the details of angular-momentum coupling, the 
single-particle level labels, and the various numerical coefficients. To highlight the re- 
arrangement taking place in each step, the factors to be reordered are indicated by an 
underbrace, and the general form of the commutator introduced by this reordering (if 
the factors do not freely commute) is shown schematically underneath. It suffices to 
note that commutation of an operator T or A through N will yield terms of the form 
[T,tf N ] _ AtStiv-i [see (42) in Scc _ 5 ] or [A,& N ] ~ S^- 1 +TS^ N ~ 1 + S^^T 
[see (43) in Sec. 5]. The time reversed adjoint expressions (see Sec. 4) apply for com- 
mutation through S N . Furthermore, new clusters will arise from commutation through 
the clusters F and G, specifically, with creation operators defined by ~ [T, G*], 
i?t ~ [T, Ft], and F ^ Q!t] ) to be calculated as described in Sec. 4. 

Beginning with the reduced matrix element = (S N G\\T\\S N F), expressed as a 
vacuum expectation value as in (15), we set out to commute the one-body operator to 



3 Relation (13) is readily derived from the identity {c\\B b \\a) = E a/ 3 7 (-) c ~ 7 ( --, a ) ( c 7\B b p \aa) [35] 
and is the operator analog of, e.g., (14.16) of Ref. [19]. We follow the normalization and phase convention 
of Refs. [30, 35] for the Wigner-Eckart theorem, i.e., {c*t\B b \aa) = (-) 2b c' 1 (aabf3\cj){c\\B b \\a) . 
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the right, where it will annihilate the vacuum. Thus, factors must be reordered as 



~ (0\(GS N ) T(S* N F*)\0) 

A t s tN-i 

~ (0\(GS N ) (S fN TIp\0) + (0\GS N A^S f N - 1 F^\0) (17) 

- (OKG^) (StJVfft)|o) + (0\GS N AW f - 1 Fl\0). 

The first term is recognized as the overlap (S N G\S N H) ~ (0\(GS N )(& N H^)\Q). Note 
that the cluster H again contains v fermions. To evaluate the second term, we must 
commute to the left, where it will eventually annihilate the vacuum, yielding 

(OlGS^AtS^-^lO) 

§N-l +T §N-l + gN-l T 

~ (0|(G^5 Ar )(^ tAr - 1 F t )|0) + (0|(G,S A, - 1 )( > S tA, - 1 F t )|0) (18) 
=1 

+ (0\{G^S N - 1 ){S^ N - 1 F^\0) + {0\(GS N - 1 )T{S^ N - 1 F^)\0) 

= E 

~ (01(75") (St^-if t)| ) + (OKGS^- 1 ) (S ,tJV - 1 .F t )|0) 

+(0|( J B5 JV - 1 )(5 tJV - 1 F t )|0) + (0|(G5 JV - 1 )T(5 tJV - 1 F t )|0). 

The first term is recognized as the overlap (S ,JV_1 (S'/)|S ,Ar_1 F), the second term 
as (S N " 1 G\S N - 1 F), the third term as (S N - 1 E\S N - 1 F), and the final term as 
(S N - 1 G\\T\\S N - 1 F). Note that each of the clusters E, H, and SI [where the S pair 
is expanded as in (7)] again contains v fermions. Therefore, the resulting recurrence 
relation has the form 

T^-O^+T^ + O^. (19) 

Similarly, a reduction in seniority v is then provided through the recurrence relation 
for the overlap 0$ = (S N G\S N F). The overlap is expressed as a vacuum expectation 
value as in (16), and a fermion is "decoupled" from each of the initial clusters F and G, 
leaving behind a subclustcr of reduced seniority, denoted by H or I, i.e., F^ ~ C^H^ 
and Gt ~ C^P. These fermion operators are then migrated inwards and recoupled to 
constitute a one-body operator, thereby yielding the reduced matrix element of a one- 
body operator, but now of lower seniority. Schematically, 

0$ ~ (0\(GS N ) (S^ N F^)\0) 

~ (0|(/G^)(^G^i/ t )|0) 

~ (0\(IS N ) (GG f ) (& N Hi)\0) (20) 
i 

~ (0\{IS N )T (& N H*)\0) + (0\{IS N ) {S^ N H^\0), 

where the latter term arises from the canonical anticommutator of C and . The result 
is recognized as consisting of a matrix element and an overlap, respectively, of seniority 
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Figure 1: Schematic recurrence network for the calculation of one-body operator reduced matrix elements 
and overlaps Ojy' in the generalized seniority scheme, as described by (19) and (21). Equation 
numbers are indicated at right for the explicit recurrence relations (Sees. 5 and 6) appropriate to each 
generalized seniority. The seed values off (circled) are given by (22). Each point (except for the seed 

values) represents multiple quantities [• • ■ ] or ojy ' [••■], distinguished by the single-particle level and 
angular-momentum coupling labels. 



v — 1. The recurrence relation thus has the form 

o^^-^ + oj;- 1 ) (2i) 



The recurrence relations of the form (19) and (21), to be set out in detail in Sees. 5 
and 6, allow any matrix element or overlap to ultimately be reexpressed in terms 
of the norms of states IS*} of generalized seniority zero, i.e., the (k = 0, 1, . . . , N). 
These constitute the seed values for the recurrence. The recurrence network for evaluating 

(v) (v) 

a given or O n is shown schematically in Fig. 1. The actual recurrence relations 
involve sums over single-particle indices and over angular momenta for intermediate 
couplings. 

The overlaps 0$ are given [3, 13], in terms of the coefficients a c defining the S pair 
in (2), by 



MeV(N.D) L 



2M, 



M 



(22) 



The sum is over all ordered partitions M = (Mi, M2, ■ ■ ■ , Md) of the integer N into D 
terms M c (J2 C M c = N), where D is the number of active single-particle levels. That is, 
the sum is over all possible ways in which N pairs may be distributed over the D levels. 



Only Pauli-allowed occupancies, M c < O c , contribute, where Q c = (2c + l)/2 is the pair 
degeneracy of the level. 

4. Commutator algebra 

In order to carry out the commutation scheme for the elementary one-body operator 

(v) 

matrix elements , as outlined in Sec. 3, it is necessary to reorder various factors 
within spherical tensor coupled products. For each value of the generalized seniority v 
being considered, it is seen from (17)— (18) that commutators of the form [A,F^] and 
[T, F-^t] will be needed, where F't is the fermionic cluster creation operator containing 
v creation operators. That is, Ft = (A^ x • • • ) x or ((A^ x • • • ) x A^) x Ct. A 
straightforward approach is provided by the coupled commutator methods of Chen et 
al. [28], as summarized in this section. All such commutators for v < 4 (some of which 
may be found in Refs. [3, 4, 22, 27, 28]) are collected in coupled form below. 
The spherical tensor coupled commutator [27] is defined by 

[A a ,B b ]t=J2(^b/3\ee)[A a a ,B b ]. (23) 

afj 

Since fermionic creation operators obey a canonical anticommutation relation, it is more 
useful in the present context to introduce the graded coupled commutator, developed by 
Chen et al. [28], in which the commutation bracket is defined by 

[A" a ,B b }=A a a B b -6 ab B b A a a , (24) 

that is, the commutator (6 a b = +1) if either angular momentum a or b is integer (quasi- 
bosonic) but the anticommutator {9 a b = —1) if both angular momenta a and b are 
odd half-integer (quasi- fermionic). [Note Q a b ^ 9{ab).) The quantity [A a ,B b ] e itself 
constitutes a spherical tensor, indeed, manifestly so, since it may be expressed in terms 
of coupled products of A a and B b as 

[A a , B b ] e = (A a x B b ) e - e ah {-) e - a - b {B b x A a ) e . (25) 

This provides the basic relation for reordering factors within coupled products. The use of 
coupled commutators thus circumvents the cumbersome process of uncoupling products 
of operators, commuting the components, and recoupling the result. The commutator 
has the symmetry (or antisymmetry) property 

[B b , A a ] e = -e ab {-) e - a - b [A a , B b ] e (26) 

under interchange of its arguments. 

The canonical anticommutation property for fermionic creation and annihilation op- 
erators, [C a , a7 Cl p] — S a bS a p (recall this is the graded commutator), is written in coupled 
form as 

[C a ,Cl] e = a5 ab 5 eQ . (27) 
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Any coupled commutator of more complicated operators can be reduced to (27) by re- 
peated application of the coupled commutator product rule 



[A a ,(C c x D d ) f } 9 = J2 



c d I 
g a h 



{{A a ,C c ] h x D d ) g 

+ ^{-) d+9 - f - h e, 



c d f 
a g h 



(C c x [A a , D d ] h ) 9 , (28) 



which applies for arbitrary spherical tensor operators A a , C c , and D d . The quantity in 
brackets is the unitary 6-j symbol, 



a b c 
d e / 



(-) 



cf 



a b c 
d e / 



(29) 



commonly denoted by U(abed; cf). 

Commutators of the one-body operator T with the creation operators for fermion 
clusters are obtained by repeated application of the product rule (28): 



6(cde)-5 bc 5 ad C\ 
d 



[T e ab , A f J} 9 = Hcd; f)9(cdegy+5 bd 

dg 



f e g 
a c d 



A9\ 



l T rsi (Kb X C t) 



9}h 



\2F(ab;e)6(abtx)^5 bs 
„ bx 



(30) 
(31) 
(32) 



e 


c 


.9 




e 


t 


X 


h 


t 


X 




r 


a 


b 



+9(ght)-6 c 



e c g 
t h r 



PL (AH x A f Jf] h = £[1 + 6(efg)nabe o cdf)]¥(cd; f) 



(33) 



x9{cdfght)^S ds 
dx 



e / 9 




7 t x 


t h x 




red 



(All x Af c )K 



The exchange symbol P(a&; e) (e.g., Refs. [3, 27]) indicates that the following expression 
should be expanded to consist of two terms related by interchange of the indices a and 

b, as 

P(a6; e) /(a, b) = /(a, b) - 6(abe)f(b, a), (34) 

and P(abc «-> xyz) indicates interchange of the indices abc with xyz. 

Commutators of the pair annihilation operator A with the creation operators for 
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fermion clusters are likewise obtained by application of the product rule (28): 



[% a ,{A<JxCt)°] h 



—8(ab; cd; e)C d 
d 

eS e fS g oS(ab; cd; e) 



-¥(cd;f)¥{ab;e)0(adef)+6 bc 

bg 

- 9(ght)^S et S ch S(ab; rs; e)C\ 
h 

+ P(rs; t)P(ab; e)9(egrt)^ts ah 6 br 5 c 
db 



e / 9 
dab 



rp9 

1 da 



(35) 
(36) 

(37) 



c e g 
a t b 



hi 



Ct+Af 



[Al s , (A e J x A U) 9 \ h = [1 + 0(efgW(abe » cdf)\U eh 5 ft 6(cd; rs; f)A n ab 

h 

-f(r S ;t)f{cd;f)f{ab;e)efgi8 br 8 ds \d \ UA^+Af, 



where we adopt the shorthand notation 

S(ab; cd; e) = 5 ac 5 bd - 0(abe)5 ad 8 bc . 



(38) 



(39) 



The 9-j symbol arises as a sum over products of 6-j symbols from (28), by identity (9.8.5) 
of Ref. [30]. In (37) and (38), Af represents additional normal-ordered terms contain- 
ing at least one annihilation operator, which are not explicitly needed for the present 
calculations, since in Sec. 5 these commutators act directly on the vacuum. 

Commutators involving the time-reversed adjoint of the creation operator for a 
fermionic cluster also arise in (18), from G 9 . However, these may be converted into 
commutators of the type considered above by use of the relation 



and therefore 



[A a ,B b f = (-) e - a - h [B h \A a ^] e . 



(40) 
(41) 



5. Recurrence relations for matrix elements of one-body operators 

The recurrence relations for are obtained by applying the commutation scheme 
outlined in Sec. 3 to the one-body operator reduced matrix element defined by (15). The 
requisite commutators are [22] 

[St»,T&t=Na d AQSt ir - 1 (42) 
[S* N ,A f cd ] f = -a^SfoNS^^-acNS^-X^ + eicdftaaNT^Si"- 1 , (43) 
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which follow from (31) and (36), respectively, by induction on N. The general recurrence 
relation 



(S N G 9 \\T* s \\S N F f ) = (-) f - t - 9 g{S N G 9 \S N H 9 ) -N 2 a 2 s (S N ^ 1 G 9 \\T t rs \\S N ' 1 F f ) 

T { N v) l--] O^l---] 

-N 2 a 2 r ff5 rs 6 fg S t0 {S N - 1 G 9 \S N - 1 Ff) i (44) 

is obtained, in terms of new fermion clusters defined by the coupled commutators 

E^ = [T^ S ,G 9 ^ 

iJft — [T* s ,F fJf ] 9 (45) 

Note that each state appearing in the matrix element and overlaps on the right hand 
side of (44) is, like the original states \S N F^) and \S N G 9 ), also of generalized seniority v. 
In particular, the cluster creation operators E* t and H 9 1 each again contain v fermion 
creation operators, and, although 1^ contains v — 2 fermion creation operators, it appears 
multiplied by a pair creation operator A®]., and the combination again carries 

(v) 

generalized seniority v. The matrix element has thus been recast in terms of an 
overlap 0$ of the same N and v, as well as a matrix element T^\ l and overlaps 0^_ 1 
of lower pair number, providing the basis for recursive calculation with respect to pair 
number. 

For v = 0: As a simple special case, the general recurrence formula (44) gives 

T#> [(aaf] = -N 2 alT£U(™) } N 2 a 2 a aO^_ v (46) 

This may readily be expressed in closed form, in terms of the seed values 0$ given 

in (22), as T^[(aa)°] = Ek=o\~) N+ HNl 2 /k\ 2 )al (N - k) aO^ . 

For v = 1: The recurrence relation obtained from (44), with the identifications F$ t — > 
C\ and G»t c], is 

TW[d\(r8)*\c] = -tS cs 6 dr O^[d\d] - N'alT^dKrsYlc] 

- N^irs^araJSdsScrO^clc] - A^rc^Ad^oO^JcIc]. (47) 

The coupled commutators needed for evaluation of P\ and H 9 ^ are given by (30) 
and (35). 
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For v = 2: The recurrence relation, obtained with — > A„\ and G 9t — > AQ , is 



[( C d)/|(r S )*|(a6) e ] 



P(afe; e)0(a&e)— £ 6s 



e t f 
r a b 



0$ [(cd)f | (ra)f] - N 2 a 2 s T£U(cd) f \ {rsf \ {abf 
^ ^ Na s a x xi5 a bSt fS e0 S(rs; cd; t)0 ( ^ ) _ 1 [(xx)°\ (aa) ] 



-N 2 a r a s ¥{cd] f)9(cdet)^S ds 
d 



f t e 
red 



O^UircYKabY 



- N 2 a 2 r re6 rs 6 ef 6 t oO { Zi[(cd) e \(abn (48) 



The necessary coupled commutators are given in (31) and (36). 

For v — 3: The recurrence relation, obtained with F^ -> (A e J x C\) 9 and G 9t 
(4] x is 



= - y"F(ab;e)6(abghx)^5 bs 
„ bx 



e i g 
h t x 



e t x 
r a b 



0^[{cdYj h \{rafi h ] 



th . 

— Oi. 

r 



e i g 
t h r 



{ p{(cdyj h \(abyr h ] N 2 alT^ 1 [(cdy j h \(r S ) t \(aby^} 
+ J2h Na s^x0{ght)hxS ft S gj S{cd; rs-J^l^xxffHabyi 3 ] 



±Na s a x P(rs;t)P(cd; f)9(fhrt)^S cg S dr 5 3S 



J f h 
c t d 



Of^xxfc^iabyiS] 



- V N 2 a r a s P(cd; f)9{cdxt)^-5 d 

„ dx 



f J h 
g t x 



f t x 
red 



of^rcTf^aby^} 



-N a r a s 9(ght)—5j S 



f J h 
t g r 



0P_ 1 [(cd) f r g \(ab) e i B ] 

- N 2 o^fg6 rB 5 gh 5 i0 OP_ 1 [(cd)fj''\(ab) e i'']. (49) 

The necessary coupled commutators are given in (32) and (37). 

For v — 4: The recurrence relation, obtained with F^ — > {A a l x A™^) 9 and G 9 ^ — > 
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(41 x A n k ?)\ is 

T^[{cdy{klY h \{r S y\(aby(i 3 rs] 



rhih , 



= + 9(emg)P(abe o ijm)]P(ij;m)9(ijm) — — £j S 

xO^ ) [(c ( i) / (fcO™' l |(a6) e (^) a; ' 1 ] 
-N'alT^l, [{cd)f (kl) n h | M* 1 [abf (tj) m 9 ] 
- \Na s a x [l + 9(fnh)P(cdf <-> kln)}hxS fg S nt 6{kl; rs; n) 



e 


m 


.9 






X 










r i 


J 



\rri9] 



y.O { y ) _ 1 [(xxf(cd)^\(abf{i 3 y 
+ Yl \ Na sU x f{rs- t)P(kl; n)P(cd; f)fnghtx6 dr 6 h 

X 

xOf^ixxficky^abYiijT 9 ] 
- V N 2 a r a s [l + 9(fnh)P(cdf o kln)]P(kl; n)6>(/Wn#fa)^< 



t g h 
I k n 
d c f 



f n h 
t g x 



n t x 
r k I 



xO^KcdyirkfoKabyiijr^} 

- N^lrgSrsSghStoO^iicdyikiy^iabyiijr^]. (50) 

The necessary coupled commutators are given in (33) and (38). 

6. Recurrence relations for overlaps 

The recurrence relations for are obtained by rearranging the factors within (16), 
as outlined in the commutation scheme of Sec. 3. It is first necessary to recouple the 
creation operators within t and t , so as to extract a single- fermion creation operator 
from each cluster, giving 

F n = (C}.xH h y f 
,* ♦ •♦ / ( 51 ) 

G n = (cixP r y. 

These equations define subclusters H h ^ and each containing v — 1 fermion creation 
operators. Then we deduce the general recurrence relation 



{S N G f | S N F*) = J2(-) f - r - h f- 1 



s r x 
h i f 



{S N P\\T? s \\S N H h } 



n(«-l)r 



+ 8 rs 5 ht (S N r\S»H> 1 ). (52) 



The derivation involves only angular momentum recoupling and the canonical anticom- 
mutator (27). Each state appearing in the matrix element and overlap on the right hand 
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side of (52) has generalized seniority v—1. The overlap 0$ has therefore been recast in 
terms of a matrix element T^' 1 ^ and overlap 0^ _1 \ both of lower seniority, providing 
the basis for recursive calculation with respect to seniority. 

For v = 1: For this simple special case, with the identifications 
G^ — > C\, the general recurrence formula (52) reduces to 



C\ and 



0$[c\c] 



(53) 



Since t and already consist of single-fermion creation operators, the "subcluster 
creation operators" H h<{ and are simply the identity operator. From the closed form 
expression for T { °\ this gives 0$[c\c] = Y,^(-) N+k (Nl 2 /k\ 2 )al (N ~ k) 0^ ) + O 



For v = 2: The recurrence relation, obtained with F^ 



A e J and G n 



(o) 

N ■ 



0^[(cd) e \{ab) e ] =^%6e)e" 1 



c a x 
b d e 



T£\d\(acy\b]+S ac 5 bd O$[b\b}. (54) 



The subclusters in (51) are H h ^ C\ and 7't -> C*. Only the t> = 2 overlaps 

0^[{aa) e \(aa) e ] (e even), [(cc)°|(aa) ], and 0$ [{ac) e \(ac) e ] or [(ca) e |(ac) e ] 
(e =/= 0) are nonvanishing, as may be shown by considering the balance of creation and 
annihilation operators in the vacuum expectation value for the overlap. Closed form 
expressions have previously been obtained [22, 23]. 

For v = 3: The recurrence relation, obtained with F^ ->■ (A e J x C\) g and G n ->• 
(A f J x is 

l(3) 



0^[(cd)'j»|(a&) e i»] 

= -E^)r 1 



e f 9 



T^[{cdy\{ijf\{abY\ +5 ij 5 ef O%>[(cdy\(aby}. (55) 



(2) r 



To decompose i* 1 -^ and G-^ according to (51), it is only necessary to reorder the coupled 
product so that the single-particle creation operator precedes the pair creation operator 
[i.e., A^xC^ -> Ctx At], which then serves as the subcluster. Thus, H h * -> (-) s_e_i ^ t 
and Pt ^ 

For w = 4: The recurrence relation, obtained with F't ->• (A^ x A^Y and G^t - 
(41 x A»t)., is 



0^[(crf)/(fcZ)"«|(a6) e (ii) TO «] 



c a x 
P Q 9 



= — ^5 1 0(abdmngq) 

pqx 

— S ac 9(bdmn) 



a b e 

m g p 

a b e 

m g p 



c d f 

n g q 

a d f 

n g p 



T^[{kl) n d%acY\{ij) m bP] 



Of[{kl) n d p \{ij) m V>]. (56) 



Here the clusters must be recoupled as A^ x A^ — > x (C^ x ) . This gives rise to sub- 
clusters H h t E P (-) p -^ m [ A ^ ] (A™ f x C\) p and ^ * EgR 9 "^" [ n ^ ] * 
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For v — 5: The recurrence relation, obtained with [(A e J x A^ f x C\\ u and 

Gf^[(AfJxA n J)»xCtr,i S 



0^l(cdy(kl) nh q u \(aby( tJ r g p u } = 



q p x 
g h u 



T { N 4) i(cdy(kir h \( P qT\(aby(ijr°} 

+ 5 pq 8 gh O ( * ) [{cd)f{kl) m °\(aby{i]) m v\. (57) 



In this case, since v is odd, the extraction of a singlc-fcrmion creation operator from 
each cluster requires only the reordering (A^ x A^) x — > x (A^ x A'), hence 
H"t _, (-J-s-P^t x A™t) 9 and r t (-)«-"-«(^t x ^tjfc 



7. Matrix elements of two-body operators 

Once overlaps and matrix elements of one-body operators have been calculated in the 
generalized seniority scheme, from the recurrence relations established in Sees. 5 and 6, 
the matrix elements of two-body operators follow in a straightforward fashion. The 
approach is based, once again, on expressing the matrix element as a vacuum expectation 
value, followed by commutation of the operators to a more convenient ordering, in which 
the various terms can be recognized as overlaps and one-body operator matrix elements. 

An arbitrary two-body operator may be decomposed as a linear combination of ele- 
mentary terms of the form {A e A x A^ cd ) w . In particular, if W w is a spherical tensor op- 
erator of angular momentum w, then the angular-momentum coupled second-quantized 
form is 

WW = I £(* + ^) 1/2 ( 1 + M^-^niWHIM'W {AH x AL) W , (58) 

abed 
ef 

where it should be noted that (A e J x A{ d ) w = -[{C\ x C\) e x (C c x C d ) f ] w , and where 
the reduced matrix element is taken with respect to normalized antisymmetrized states 
|(a6)g)NAS = (l + ^o6) _1 ^ 2 (C] x C^)||0). We note this explicitly to avoid ambiguity with 
the unnormalized N = 0, v = 2 basis state of Sec. 2. The expression for the two-body 
Hamiltonian arises as a special case, with w — 0. 

Evaluating the matrix element of an elementary two-body operator {A e a l x A^ d ) w as a 
vacuum expectation value involves only a single commutation, of the form [A, A^] ~ T+l 
[see (36)]. Schematically, 

{S N G\\A*A\\S N F) ~ (0\(GS N ) (A^A) (S* N F*)\0) 

T+l 

~ (0\{GS^A) (AtS^Ft)\0) + (0\(GS N ) T (S* N F*)\0) + (0\(GS N ) (S* N Ft)\0). (59) 

That is, the pair creation and annihilation operators A^ and A, which together constitute 
the two-body operator, must be decoupled from each other and reassociated with the 
fcrmion cluster operators for the two states. If F and G are clusters of generalized 
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seniority v, then the first term above is recognized as an overlap involving two new 
clusters, of higher generalized seniority v + 2. Let us represent the reduced matrix 
elements of the fundamental two-body operators between generalized seniority basis states 
of equal generalized seniority by • •], e.g., 

W ( N 4) [(cd) f (kl) nh | (rsf (xyf w | (abf {ij) m «] 

ee (S N (cd) f (kl) nh \\{A t r l x Al v T\\S N {aby{ijT°). (60) 



Then the resulting expression for the two-body operator matrix element is of the form 

IfW-O^+TW + oW. (61) 
The full relation obtained following this commutation scheme is 



{S N G°\\(A% x Al y r\\S N Ff) = j2(~y +f - k k 



z t w 
f 9 k 



(S N I k | S N H k ) 



->(" + 2)r 



tz 



+ ¥(xy; z)F(rs; t)6(sxw) — 6., 



rw 



ry 



t Z W 

x s r 



(S N G 9 \\T™ x \\S N F f ) 



- fiS tz S fg S wQ S(rs;xy;t) (S N G^\S N F^, (62) 



■>Wr 



in terms of new clusters defined by 



ff fet = (AH x Fn ) k 
I kf = (A z x l x G9t ) fe - 



(63) 



To obtain a useful relation, in terms of known overlaps of seniority v + 2, it may be neces- 
sary to recouple the factors making up the operators H k t and I k \ so that these fermion 
cluster creation operators are of the form used in defining the generalized seniority states 
in Sec. 2. For instance, for the two-body operator matrix elements involving states of 
v = 3, (S N I k \S N H k ) matches the definition of Ojp in (5) only after the recoupling 
[At x (At x Ct)] ->• [(At x At) x Ct]. 

For reference, let us explicitly write the relations for two-body operator matrix el- 
ements for states with generalized seniority v < 3. These expressions involve the T N 
with v < 4 and 0$ with v < 5, as considered explicitly in Sees. 5 and 6. 

For v = 0: 



(v) 



< 0) [(r S )W 



iO^[(xyY\(r S f] 



F(xy;t)¥(r S] t)U rx 5 sy T^[(rr) } 



i6(rs;xy;t)0fp . 



(64) 
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For v = 1: 



= -^6{ctk)k 



z t w 
c d k 



< £ ) [(xyyd k \(rs) t c k i 



tz 



+F(xy; z)¥(rs; t)6(sxw)j^5 ry 



t z w 
x s r 



T^[d\(sxr\c] 



ciS tz 6 cd 6 w0 5(rs; xy; t)0 ( ^ [c\c]. (65) 



For v = 2: 
W^[{cdy\(r S y{xyf w \{abY 



^9{etk)k 



z t w 
e f k 



0f[{xyY{cdy k \{r S f{ab) 



tz 



+¥(xy; z)V{rs; t)6(sxw) — S ry 



t z w 
x s r 



T^[(cdy\(sxr\(abY 



- ei5 tz 5 ef 5 w0 5(rs;xy;t)O^[(cd) e \{ab) e }. (66) 



For v = 3: 



[{cd) f j h \ (rsf (xy) z w \ [abfi 9 ] 



Y^0{tgk)k 

pqk 



z f q 

j k h 



t e p 
i k g 



z t w 
g h k 



OffKxyyicdyfKrsfiabyvi 



tz 



+¥(xy; z)P{rs; t)9(sxw) — S ry 
rw 



t z w 
x s r 



TW[(cdyj h \(sxr\(abyiS] 
- gtS tz S gh 6 w0 6(r S ; xy; t)0^ [(cd) f f | (ab) e i g ] . (67) 



8. Conclusion 

The calculational framework presented here provides a straightforward and system- 
atic approach to constructing recurrence relations for matrix elements and overlaps in a 
generalized seniority scheme. Matrix elements of the one-body and two-body multipole 
operators have been considered explicitly here, but matrix elements of other operators 
of interest, such as the pair transfer operator [4, 18], may be derived similarly. 

Aside from some elementary angular-momentum recoupling, the derivation only re- 
quires calculation of commutators of the form [T, x • • • x A<] (or [T,A^x---xA^x Ct] 
for odd particle number) and [A, A 1 x • • • x A 1 } (or [A, A* x • • • x A^ x C t ] for odd particle 
number). This may be accomplished systematically via the coupled commutator product 
rule (28). Although the process becomes increasingly laborious for larger v (i.e., v = 5, 
6, and 7, or higher if needed), it is also well-suited to automation [29]. 
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